Optical study of superconducting Ga-rich layers in silicon 
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We performed phase-sensitive terahertz (0.12 - 1.2 THz) transmission measurements of Ga- 
enriched layers in silicon. Below the superconducting transition, y™ lddl ° = 6.7 K, we find clear 
signatures of the formation of a superconducting condensate and of the opening of an energy 
gap in the optical spectra. The London penetration depth, A(T), and the condensate density, 
n B = A 2 (0)/A 2 (T), as functions of temperature demonstrate behavior, typical for conventional su- 
perconductors with A(0) = 1.8 /im. The terahertz spectra can be well described within the frame- 
work of Eliashberg theory with strong electron-phonon coupling: the zero-temperature energy gap 
is 2A(0) = 2.64 meV and 2A(0)/fc s T c = 4.6 ± 0.1, consistent with the amorphous state of Ga. At 
temperatures just above T c , the optical spectra demonstrate Drude behavior. 

PACS numbers: 74.25.Gz, 74.70.Ad, 74.78.-w, 74.25.nd 
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II. EXPERIMENT 
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Recently, some of us observed superconductivity in 
amorphous Ga-rich layers manufactured by Ga implan- 
tation into Si wafers and subsequent thermal annealing X 
Compared to bulk crystalline Ga, superconductivity oc- 
curs at higher temperatures with an onset at 7 - 10 K. 
While previously, amorphous Ga films have only been 
characterized subsequently to in situ low-temperature 
synthesis the films under investigation withstand mul- 
tiple cooling procedures and room-temperature handling. 

Here, we report optical evidence for the 
supcrconducting-statc formation in gallium-enriched 
layers in silicon by means of optical measurements in 
the terahertz region. These layers consist of amorphous, 
gallium-rich precipitates embedded in nanocrystallinc 
silicon^ Optical spectra of superconductors contain 
valuable information about the superconducting state: 
the London penetration depth, the strength of coupling, 
the size and the symmetry of the superconducting gap 
can be extracted from such measurements 

In our optical measurements within the terahertz re- 
gion, wc observe clear signatures of the superconduct- 
ing condensate developing at T < T c . Wc were able to 
trace the temperature dependence of the spectral weight 
of the condensate and of the London penetration depth. 
We demonstrate that they follow nicely the behavior 
expected for fully gaped superconductors. We further 
show that the frequency-dependent optical spectra can 
be well described within the Eliashberg theory for strong- 
coupling s-wave superconductors. 



Ga + ions had been implanted at an energy of 80 keV 
with a total fluence of 4 x 10 16 cm~ 2 Ga into 0.38 mm 
thick (lOO)-oriented silicon wafers covered with 30 nm of 
silicon dioxide. Subsequent rapid thermal annealing at 
650 °C for 60 seconds has been applied for realizing gal- 
lium precipitation at the Si - Si02 interface. The thick- 
ness of the fabricated Ga-rich layers was estimated to 
be 10 nm by means of transmission electron microscopy 
and Rutherford backscatteringi^ The samples probed op- 
tically had lateral sizes of 2 by 5 mm. Direct-current 
resistivity measurements (Fig. [1} show the onset of the 
superconducting transition at T° nsct « 7.3 K and a mid- 
dle point at T c middlc = 6.7 K, which we take as T c for our 
Eliashberg analysis below. Further information, such as 
a detailed structural and critical-field analysis, can be 
found in Rcfs. Hand!. 

In the frequency range 3.9 - 41 cm" 1 (117 - 1230 GHz, 
0.48 - 5.1 meV), measurements have been performed by 
use of a spectrometer, equipped with backward-wave os- 
cillators (BWOs) as sources of coherent and frequency- 
tunable radiation. The measurements have been done 
using a number of different BWOs covering the above- 
mentioned frequency range almost continuously. 

A Mach-Zchnder interferometer arrangement of the 
spectrometer allows to measure both the intensity and 
the phase shift of the wave transmitted through the sam- 
ple. A commercial optical 4 He cryostat with the sample 
inside is installed in one of the arms of the interferome- 
ter. A mirror in another (reference) arm is connected to 
a precision stepper motor (with an accuracy of 0.5 /im). 
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FIG. 1: (Color online) Temperature dependence of the resis- 
tance of a Ga-enriched layer in Si (right-hand scale) and an 
example of the raw measurements of the interferometer mirror 
position (left-hand scale). The mirror position is directly re- 
lated to the optical phase shift and, hence, to the permittivity 
and the penetration depth. 



The interferometer can be adjusted to a position, where 
the optical path difference between the arms is zero. Dur- 
ing frequency or temperature sweeps, a software always 
keeps the movable mirror in position with zero optical 
path difference. Any change in the sample's refraction in- 
dex (or dielectric permittivity) results in a change of the 
optical path length in the sample arm, and, hence, to a 
shift of the movable mirror. This shift is detected and the 
phase change of the wave transmitted through the sample 
is directly calculated from it. For the amplitude trans- 
mission measurements, the reference arm of the interfer- 
ometer is blocked with a shutter. The absolute values 
of the sample's amplitude and phase-shift transmission, 
Tr and tp, are obtained by repeating the measurements 
without a sample and a subsequent correction for the 
empty-channel measurement. Fresnel optical formulas^ 
are used to extract the optical parameters of the sample 
(for example, the complex conductivity, a = <7i + i<?2, 
or the complex permittivity e = e' + is") from Tr and 
tp. A detailed description of the measurement technique 
can be found in Ref. |8j. The method has been previously 
applied to a large number of different superconductors^ 

The London penetration depth, A = cj {^a^tS) 1 ! 2 {to 
is the angular frequency, c is the speed of light), has been 
determined as a function of temperature at a number of 
fixed frequencies (123, 172, 217, and 252 GHz, or 4.1, 
5.74, 7.24, and 8.4 cm -1 ) in the same way as described 
in Refs. [lfl and [nj For these measurements, we limited 
ourselves to our lowest frequencies, as at higher frequen- 
cies the contribution of normal electrons to a 2 becomes 
significant, thus, a correct determination of A is difficult. 

The optical parameters of pristine wafers had been 



measured in advance using the same setup. At T < 15 K, 
we found the real part of the dielectric permittivity of the 
wafers to be independent of temperature and frequency 
(within our accuracy) and equal to e' = 11.50 ± 0.05. 
An absorption in the wafers was not detectable at these 
temperatures. 



III. RESULTS AND DISCUSSION 

A. Temperature sweeps. Temperature dependence 
of penetration depth and condensate density 

Firstly, let us note that the appearance of superconduc- 
tivity in our sample is confirmed by the raw optical data. 
In Fig. [TJ we show an example of such raw data - the 
position of the movable mirror of the Mach-Zehnder in- 
terferometer, x(T), as a function of temperature. Above 
the superconducting transition, x(T) is flat. As the sam- 
ple enters into the superconducting state, x(T) starts to 
decrease rapidly. This is due to the condensation of elec- 
trons. The condensed electrons are represented by a delta 
function in the real part of the complex conductivity, 
01. This delta function leads (via the Kramcrs-Kronig 
relations) to a divergence in 02 (and in e') at to — > 0: 
a 2 = ne 2 /mjjj, e 1 = 1 — Air^/u = 1 — 4:irne 2 /mui 2 (n 
is the charge-carrier concentration, e is the elementary 
charge, and m is the effective mass of the carriers). One 
can easily estimate that for conventional superconduc- 
tors at temperatures well below T c , the permittivity at 
(sub)terahcrtz frequencies is of the order of — 10 3 to — 10 6 . 
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FIG. 2: (Color online) London penetration depth as a func- 
tion of temperature for some frequencies. Inset: Relative den- 
sity of superconducting condensate, n s = A 2 (0)/A 2 (T), as a 
function of temperature for the same frequencies as in the 
main panel. Solid line mimics the behavior expected for a 
fully gapped superconductor. 



3 



This negative e' leads to a negative contribution to the 
phase shift and hence to a negative contribution to x(T), 
directly recorded in our measurements. Obviously, this 
huge negative e' is nothing but screening of the prob- 
ing electromagnetic field by supercurrents. As electrons 
start to condense already at T° nset , x(T) begins to de- 
crease starting at this temperature. 

The measured temperature dependences of A for the 
above-mentioned frequencies are shown in Fig. [5J All 
curves demonstrate a very similar behavior. At low tem- 
peratures, the curves flatten. Extrapolation of A(T) to 
zero temperature gives A(0) = (1.8 ± 0.2) /mi. As the 
temperature approaches T c , A(T) nearly diverges for all 
four frequencies (because at any non-zero frequency, <72 
is different from zero above T c , true divergence is not 
expected). 

The inset of Fig. [2] shows the relative density of 
the superconducting condensate, n s = A 2 (0)/A 2 (T), as 
a function of temperature. All measurements basically 
fall onto one curve, which can be well approximated by 
a standard phenomenological expression used to fit the 
penetration depth data of conventional (fully gapped) 
superconductors £ n s = 1 — (T/T c ) 4 . The highest- 
frequency curve deviates slightly from the common trend 
as T — > T c . We attribute this deviation to the contri- 
bution of normal-state carriers, as mentioned above (a 
similar behavior was observed e.g. in Rcf. 10). Our 
measurements are more noisy at lower frequencies be- 
cause the wavelength of the probing radiation at lower 
frequencies (123 and 172 GHz) becomes comparable to 
the lateral sample dimensions. 



B. Frequency sweeps. Drude metal and 
strong-coupling s-wave superconductor 

Normal state. Figure [3] shows (a) the as-measured 
transmission, Tr(i/), and (b) the phase-shift, f{v), spec- 
tra as a function of frequency, v = lo/2tt at 10 K. Since 
the major term of the phase shift is proportional to the 
frequency of the probing radiation, the phase-shift spec- 
tra are divided by frequency to eliminate the constant 
frequency slope. The pronounced fringes in both Tr(i/) 
and tp(v) are due to the multiple interference within the 
Si substrate, which acts as a Fabry-Perot interferometer. 

The spectra can be well fitted by the Drude model, 
cr{tjj) = Co/(l ~* wt )j where Co is the zero-frequency limit 
of the conductivity and r is the scattering time. We find 
(To = (900 ± 150) f2 _1 cm _1 in good agreement with our 
dc data [cr dc (10K) = 1000 n _1 cm -1 ]. The scattering 
rate, 7 = 1/(2ttt), is above our frequency range. We 
can estimate, 7 > 100 cm -1 (r < 5.3 x 10 -14 s). Lower 
values of 7 would result in an upturn in Tt(v) at the high- 
frequency end of our range, which we do not observe here. 
By use of the Fermi velocity for gallium, vf = 6 x 10 7 
cm/sj 1 ^ we can estimate the upper limit of the mean free 
path at 10 K, I < 30 nm. 

At our high frequencies, we do not see a downturn 
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FIG. 3: (Color online) Frequency-dependent amplitude- 
transmission (a) and phase-shift (b) spectra of the Ga-rich 
layer in Si at T = 10 K. Solid lines show a Drude fit. The 
phase-shift spectra are divided by frequency for better repre- 
sentation. 



of transmission either. Such a downturn would indicate 
greater o\ values at higher frequencies and, hence, signal 
spatial localization of the carriers. Instead, our Drude- 
type conductivity proves a free-electron transport in the 
Ga layer. From the value of our highest measurement 
frequency, we can roughly estimate that no localization 
happens on scales of a few dozens of nanometers or larger. 
One cannot, of course, exclude localization effects on 
smaller spatial scales. 

Superconducting state. The changes in the spectra re- 
lated to superconductivity can best be seen in comparison 
with the normal-state spectra. Therefore, in Fig. 21 the 
ratio between the transmission coefficients at 4 and 10 K, 
Trg/Trjv, is shown in panel (a); the difference between 
the phase shifts, ips — ^Pn, measured at the same temper- 
atures, is in panel (b). The presence of fringes in both 
panels indicates a large change of the dielectric constant, 
which occurs when the sample enters the superconduct- 
ing state. 

In an attempt to compare experiment with theory we 
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FIG. 4: (Color online) Transmission ratio, Trs/Trjv, (a) and 
relative phase shift, tps — <Pn, (b) for the Ga-rich layer in Si 
as functions of frequency. Normal-state data are from Fig. [3] 
Superconducting-state data are collected at 4 K. Solid lines 
show the Eliashberg calculations for T = 4K with 2A(4K) = 
2.55 meV = 20.6 cm -1 . The phase-shift spectra are divided 
by frequency. The real, <Ji(v), and imaginary, 02 (f), parts of 
the complex conductivity obtained from the calculations are 
shown in the insets. 

calculated the complex conductivity within the frame- 
work of standard s-wave Eliashberg theory ! 13 i 14 Details 
of the procedure applied in these calculations will be dis- 
cussed elsewhere In agreement with our data on the 
temperature dependence of the penetration depth and 
with the amorphous state of Ga in our samples^ we as- 
sumed an isotropic s-wave gap. The Eliashberg function, 
a 2 {uj)F{uj), and the Coulo mb p seudopotential, fi* = 0.17 
have been taken from Ref. [LQ. The Eliashberg function 



was rescaled to give T™ lddlc keeping fi* unchanged. This 
results in a mass-enhancement factor due to electron- 
phonon interaction of X e i- P h = 1-97 down from 2.3 as 
was reported by Chen et alJ£ We obtained the zero- 
temperature energy gap: 2A(0) = 2.64 meV = 21.3 
cm -1 , and 2A(0)/k B T c = 4.6 ± 0.1 (the error is esti- 
mated from the width of the superconducting transition). 
This value of 2A(0)//cbT c is in good agreement with the 
results of earlier tunneling and infrared experiments on 
amorphous Ga»iZr— 

The calculated complex conductivity at T = 4 K is 
shown in the insets of Fig. H](b). The theoretical complex 
conductivity is then used to calculate the transmission 
ratios and the phase-shift difference via Fresnel formulas 
[lines in Fig. [4]. The good agreement between theory 
and experiment is evident. 

Finally, from the obtained value for 2A(0), we can 
estimate the zero-temperature coherence length: £0 = 
?ivf/ttA(0) = 100 nm. Let us note that although the 
Ginzburg-Landau coherence length £* derived from the 
analysis of the upper-critical-field data, was found to be 
around 10 nmji there is no contradiction with our result, 
because £*, derived in that way, is strongly affected by 
the mean free path. 



IV. CONCLUSIONS 

Our optical data collected at terahertz frequencies on 
Ga-enriched layers in Si: i) support the original ob- 
servation of a superconducting transition in resistivity 
data 1 - via the direct detection of the electromagnetic- 
field screening due to supercurrents; ii) show that the 
London penetration depth [A(0) = 1.8 /im] and the su- 
perconducting condensate density as functions of tem- 
perature follow the typical temperature dependence ex- 
pected for conventional superconductors; iii) prove that 
the frequency-dependent normal-state conductivity is of 
Drude type (free electrons) with no signs of localization 
effects on length scales equal or larger than several dozens 
of nanometers; iv) indicate that the upper limit of the 
mean free path £ in our sample is roughly 30 nm; v) allow 
to estimate the size of the superconducting energy gap, 
2A(0) = 2.64 meV, and the coherence length, £ = 100 
nm, and to conclude that our sample is a dirty (£0 > 
strong-coupling (\ e i- P h = 1-97) s-wave superconductor 
with 2A(0)/fc B T c = 4.6 ± 0.1. 

Part of this work was supported by EuroMagNET II 
(contract 228043) and by the DFG (contract HE 2604/7- 
1). 
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